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An exact solution was recently found in the massive gravity theory having the form of 
Schwarzschild-dS black holes with some additional background fields. Hawking radiation will occur 
at the event and cosmological horizons having the blackbody spectrum, which will be modified by 
the geometry outside the black hole. In this paper, we study the greybody factors of a test scalar, 
considering its minimal coupling with the background geometry. The case of small black holes with 
the horizon radius much smaller than the cosmological dS radius is studied numerically. The case 
of near-extremal black holes with the horizon radius comparable to the cosmological dS radius is 
studied analytically. In addition, we considered the coupling of the test field with the background 
Stiickelberg helds, which in turn leads to reductions in particle emission and some non-trivial features 
(resonances) in the greybody factors. 

PACS numbers: 


I. INTRODUCTION 


The Lagrangian we will consider is 


Hawking emission from a black hole is characterized 
by two factors: the black body factor which gives the 
probability that a particle is created in the vicinity of a 
horizon, and the greybody factor which gives the prob¬ 
ability that this particle penetrates the potential barrier 
and escapes to infinity. The pedagogical difference be¬ 
tween these two factors can be found in [1] (see also 0). 
The number of distinct models with the exact black hole 
solutions is constantly increasing, which in turn implies 
the need for concrete calculations of the corresponding 
greybody factors (for some recent developments see for 
example [SHU]). The main goal in this paper is to cal¬ 
culate the greybody factors for black holes in massive 
gravity, which is the case that has not been studied in 
the literature so far. 

Massive gravity is a collective name for theories where 
gravitons appear to be massive. The original idea is old, 
appearing for the first time in m. but gained a lot of at¬ 
tention recently. Namely, massive gravitons emerge nat¬ 
urally in some higher-dimensional models like the Dvali, 
Gabadadze, Porrati model cascading gravity [181- 
124] and ghost-free massive gravity m- For nice reviews 
please see [13 US] . 

In massive gravity models, the asymptotically flat so¬ 
lution of vacuum Einstein equations converts the horizon 
to a spacetime singularity. One way out of this problem is 
to study the asymptotically de Sitter solutions |2H| • Then 
the simplest case is a Schwarzschild-dS black hole, with 
background Stiickelberg scalar helds. As a Hrst study of 
Hawking emission in this class of black holes, we consider 
the emitted scalar held minimally coupled to the gravi¬ 
tational background with and without coupling to these 
external helds. Though emission from a Schwarzschild- 
dS black hole has been studied before, coupling to the 
background helds will introduce some non-trivial features 
in the greybody factors. 


/ 1 \ 

C = + d^Lpd,,ip\ . 

( 1 ) 

Here rrig is the graviton mass and lA is the graviton 
potential which depends on the metric as well as the 
Stiickelberg helds 4’°'. Explicitly, 

U{g, r) =U, + -(^- + U,, (2) 

where a^, 04 are free parameters, and 
U2 = 

U 3 = [/C]3-3[A][A2] + 2[/C3], 

- 6 [A 2 ][A]^-b8[A3][A]-b3[A2]2- 6 [A‘^], 

(3) 

with the matrix K,^[g, 4°') = (5^ — \/ gi^°'da4°'d,y4’’Vab, and 
[/C] is the trace of it. 

As found in [55] , neglecting the backreaction of the test 
helds If on spacetime, there exist Schwarzschild-dS type 
of black hole solutions for a special choice of parameters 
a and 77 as ry = — ^. 


ds^ = -(I— --i^rug 
r 6a 


^)df + 


dr^ 


I- - — 


.+r^dn^ 


3cx 


, (4) 

where the free parameter Vg is the gravitational radius of 
the black hole. In this solution, the cosmological horizon 
size is determined by the the graviton mass. In what 
follows, we consider the equation of motion of ip in this 
given spacetime. 
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II. GENERAL ANALYSIS 


Now we consider the spherically symmetric spacetime, 
ds^ = -F{r)dt^ + F-\r)dr^ + r^dVt^. (5) 


For the scalar Hawking emission, the equation governing 
its kinematics is the Klein-Gordon equation. The scalar 
field must retain the form where the variables can be 
separated as 

( 6 ) 

r 


where w is the (real) frequency under consideration, and 
Yim are the spherical harmonics. Now the Klein-Gordon 
equation reduces to a single radial equation. 


' d^ 
dr*'^ 


Vir) 


= 0 . 


(7) 


for near-extremal black holes Rh ~ Rc- 

We can see that F{r) vanishes at both event and cos¬ 
mological horizons, and that r* —>■ — oo as r —>■ RH^ 
r* —>• - 1-00 as r —>• Rc- Therefore, the solution to Eq. Q 
is readily written once the boundary condition is speci¬ 
fied. 

First, if particles are emitted from the event horizon, 
travel to the cosmological horizon and are partly reflected 
back, then the solution (”up” modes) assumes the follow¬ 
ing form. 


^ -f ^-oo 

~ ^-koo. (16) 

Second, if particles are injected from the cosmological 
horizon, partly reach the event horizon and partly get re¬ 
flected back, then the solution (”in” modes) should have 
the form 


Here r* is the tortoise coordinate defined as 

* - f 

J F{ry 


( 8 ) 


and V(r) is the effective potential determined by the met- 
ric, 


V(r) 


F{r)drF{r) 

r 


+ l{l + l) 


F{r) 

J.2 * 


(9) 


^ ^-oo 

~ -k , r* ^ -koo. (17) 

For real frequencies, the radial equation is invariant un¬ 
der complex conjugation, so there is another solution, 

^-oo 

^ ^+oo. (18) 


We now rewrite F{r) concretely for Schwarzschild-dS 
black holes in the popular form, 

F{r) = 1 — Tg/r — (10) 


with K = \l ^nig. The event horizon Rh and the cos¬ 


mological horizon Rq will be two of the three roots of 
F{r), and the third one denoted as Rx, i.e. 


Fir) = 


K^ir - i?_f/)(r - Rc)ir - Rx) 


( 11 ) 


By extracting out the constant term in the numerator, 
we get 


Tg = —RhRcRx > 0 , ( 12 ) 

thus Rx = — Hh^r ^ 9, which means that there are no 
singularities, coordinate or physical, between the event 
and cosmological horizons, which makes it straightfor¬ 
ward to numerically solve Eq. Q between Rh and Rc- 
In terms of the parameters in massive gravity. 


Rh 


Rc 


v/3^ 

mg 


for small black holes Rh «C Rc, and 


Rh ~ Rc 


_ 3rg 

mg 2 


(13) 

(14) 


(15) 


A linear combination of solutions ( [T 6 | and (18 1 is also a 
solution of the linear equation Q, so we can define 


4- - (l-i?i?*)e*‘^’' ,r* ^ -oo 


-koo. 


(19) 


Equations 0 and ( [T^ are two linearly dependent solu¬ 
tions, thus 


T' 1 R' 

1 - RR* ~ jF - ~ ji*T’ '' ' 

Combining this with the definition of the greybody fac¬ 
tors 7 ( 0 ;) = TT* = 1 — RR*, we get 

This means that the black hole particle transmission co¬ 
efficient (greybody factor) is equal to the transmission 
coefficient for emissions from the cosmological horizon. 

In Penrose diagram. Fig. Q |53], particles emitted 
from the past event horizon 'H~ partly travel to the fu¬ 
ture cosmological horizon Hq and partly get reflected 
back to the future event horizon 77^. Similarly, particles 
emitted from the past cosmological horizon Ry are partly 
received at the future event horizon , and partly re¬ 
flected to the future cosmological horizon Hq. And as 
shown above, the greybody factors for these two pro¬ 
cesses are equal. 
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FIG. 1: Partial Penrose diagram for the 
Schwarzschild-dS spacetime, along with the emission 
and transmission of the emitted particles from both 
horizons. [29] 


III. NUMERICAL RESULTS 


In this section we restrict our attention to small black 
holes, that is, Rh «C Rc or M <C 1 /k. We chose kRh = 
0.01 without loss of generality. For this horizon radius, 
the surface temperature Th = ^ = ^ ~ 8k. To 

obtain the black hole emission rate and spectrum, we 
take the cutoff I to be 3. We adopt the fourth-order 
Runge-Kutta method in the numerical integration. To 
get convergent results independent of the computation, 
we chose the step size of integration as 10“^ for Z = 0,1, 2, 
and 10“"^ for Z = 3 to treat the higher frequency modes 
correctly. The results are shown in the figures (1|)(§(14|. 

As seen from equation the potential barrier gets 
higher as I gets larger, making it more difficult for 
low-energy particles to get across, consistent with Fig. 
(H). Besides, the contribution to the total emission rate 
and power from modes with frequencies well beyong the 
Hawking temperature is exponentially depressed. So we 
restrict ourselves to uj/k < 100, at which energy the mode 
I = 3 contributes less than 1/1000, as seen from Fig. 

This gives roughly a precision of one part in a thousand 
in the spectral emission rate and power spectrum we ob¬ 
tained in figures 

Our result for Z = 0 at low frequencies is consistent 
with the analytical approximation of Harmark et al. |29j . 

7oM«4(Ki?ff)2(l+ccVAt"), (22) 


as was derived by matching the solution forms in different 
regions of r, and is also consistent with the result of grey- 
body factors in the w —> 0 limit for (4 -|- n)-dimensional 
Schwarzschild-dS black holes by Kanti et al. [55] , 


70 ,n(w — 0) — 4 




(23) 


in the case of n = 0 and Rh ^ Rc- Qualitatively, our 


numerical results also agree with those by Crispino et al. 
[35] . though for a different small event horizon size. 

We also plot the expected emission rate {n{uj)) and the 
power spectrum , which have the well-known expres¬ 
sions as 


(n(a;)) 

dE 

dtdX 


E 


(2Z -I- l)7z(a;) 
exp(27ra;/fc//) — 1 


(24) 


i 

dE 
2 tt dtduj 


E 


(2Z -I- 1 ) 7 ; (w) 


^ exp{2TTUj/kH) — 1 


(25) 


Compared with the blackbody radiation, the emission 
rate is largely reduced at low frequencies due to the small 
greybody factors there. This reduction shrinks at higher 
frequencies. Also, the peak frequency on the power spec¬ 
trum moves to higher values, from that of the Planck’s 
radiation specturm. 



FIG. 2: Greybody factors for Schwarzschild-dS black 
holes for Z from 0 to 3. The thin vertical black line 
shows the position of the Hawking temperature, and the 
thick solid black line is the low-frequency appximation 
in Eq. (22). 


IV. NEAR-EXTREMAL BLACK HOLES 


In this section we consider the case of near-extremal 
Schwarzschild-dS black holes, where the black hole event 
horizon and cosmological horizon are very close. Now 
we can get approximately a closed form of V{x), to the 
leading order in surface gravity kn « Then, 


1 


This is the Poschl-Teller potential, which makes the 
Schrodinger equation exactly solvable [3D|. With the 
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FIG. 3: The expected emission rate vs. frequency, 
including modes with / = 0,1, 2,3, for small 
Schwarzschild-dS black holes. 



FIG. 4: Emission power spectrum of small 
Schwarzschild-dS black holes, compared with the 
blackbody power spectrum. Here A = 27r/w, where we 
have taken the speed of light to be I. The greybody 
spectrum is amplified by a factor of 10 to be more 
visible. dE/dtdX has the unit Note that the 
greybody spectrum is always below the blackbody one 
before the artificial amplification. 


where R and C can be written in terms of Gamma func¬ 
tions [30]. The greybody factor is readily obtained as 

^ |sin(7r/3)r(-/3 + *^)r(l + ^ + z^)|2 

i 2 • 


Note that this result doesn’t hold for ^ = 0, in which 
case the potential V{r) ^ thus the transmis¬ 

sion coefficient will be much higher than the I > 0 cases. 
On the other hand, for very large I and low frequencies 
around the scale of kn, 7i(w) vanishes, from the property 
of gamma functions. 


n-p+i—ni- 

kn 


■ ^ \ 


r(-/3)F(l+/3) 


sin(7r/3) ’ 


(30) 


At high energy and I ^ 10, then the arguments in all 
functions in Eq. (29) are practionally imaginary, mak¬ 
ing further simplification possible. Using the asymptotic 
approximation 


27r 


1/2 


r(/x)«e-“(zx)“ — (l + 0(l/a;)), (31) 

' IX ' 


we can put 7 /(w) in a cleaner form for w not close to 
, so that arguments of the gamma functions in 


Eq. ( |29| ) are both imaginarily large. Therefore, keeping 
only the leading order term. 


7i(w) 




1 - 


/ tjJ I I \( ^ _ / \ 

•V /^t:- 7 ? XX ' ' fc •V rr ' 


= 1- 


4e 


\/^kRh \/3k.Rh 

— 27r max 


■ sinh 


2 f a 

V V^kRh / 




(ft)' 


(2 




sinh^ 


V VskRh J 


Thus for w > , 

VSkRh 


(32) 


7 i(a;) « 1 - ( — j e 


^h\ ^-2TTUj/kH 


(33) 


boundary condition ik ^ as r* —>■ — oo, we can 

get 

'I'(^*) = K(1-C)]^2Ai(1 + ^+^,^-/3;1 + ^;C), 

kn kn kn 

(27) 

with c = (1 +p = _ i_ As 

r* —>• - 1-00 or ^ I, the above solution has the asymp¬ 
totic form 

4'(r*) ~ C (e”’''' -f , (28) 


which is exponentially approaching I as w increases. And 


7z(w) 


HkRh)^ 

P 


(34) 


which is independent of w, but approaching I quadrati- 
cally as I increases. 

For I from I to 8, our analytical expression will be 
shown to be perfectly consistent with the numerical re¬ 
sults in figure 
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V. COUPLING WITH BACKGROUND FIELDS 


As we have demonstrated in previous sections, without 
the coupling of the test scalar field to the background 
Stiickelberg fields, calculations effectively reduce to the 
Schwarzschild-dS black hole case. In this section, we take 
the coupling into consideration. Massive gravitons can 
be decomposed into tensor modes, as well as vector and 
scalar modes and tt, which turn out to be the non¬ 
unitary parts of the background Stiickelberg fields [25] . 
i.e. — (j)^ = {mgA^ -\- 9^7r)/A^. In the Schwarzschild 
coordinates, the vector modes were found to be 

A* = 0, (35) 

A° = -^^/(r), (36) 

Ko 

and the scalar mode tt depends on both r and t. Here 
the free dimensionless parameter kq is an integration con¬ 
stant when solving the Einstein equation to obtain the 
Schwarzschild-dS metric [25], and 


fir) = ± 


dr 



2 9 2 

3a 


9 o o ' 
3a 


(37) 


which diverges at both horizons. 

The TT mode depends on the coordinate time [25j, so 
coupling with it will destroy our assumed solution form 
in Q. So for simplicity, we focus on the interaction with 
the vector modes A^, and further require invariance un¬ 
der the Lorentz and discrete ip ^ —Lp transformations, 
and renormalizability. This leaves only y/—gA^p^ and 
y/—gpdg,pA>^ couplings in the lagrangian density. How¬ 
ever, 


(38) 

where the first term on the right-hand side is a total 
derivative thus not contributing to the equations of mo¬ 
tion of p, and the second term vanishes because and 
the only nonzero component of the vector mode A*^ do 
not depend on the coordinate time. Therefore, from now 
on, we consider only the interaction term 

^coupling — ^V^(AAV"). (39) 


Here A is a dimensionless coupling constant. This form 
of coupling is the only possible renormalizable form con¬ 
sistent with the symmetries in the model. It would break 
any eventual gauge invariance for the field A^, since it 
will act as an effective mass for the field A^, but in mas¬ 
sive gravity this field is massive anyway. 

This coupling term induces an additional potential 
term in Eq. 


^coupling ) 


—AA^(r)F(r) 




F{r) 


(40) 


dr 


\/l - F{r) 

F{r) 


The integral term has logarithmic divergences at both 
horizons. However it is multiplied by the F{r) term in 
front of it which fixes this divergence, thus enabling us to 
use the techniques in section H to compute the greybody 
factors in this case. The forms of the effective potentials 
are shown for the small black holes in figure ([^. Here, the 
contributions from black hole geometry, i.e. Eq. © , are 
seen as the spikes on the left side of the figure, while the 
contributions from the coupling with massive gravitons 
are shown as the spikes on the right part. 

The total effective potential for near-extremal black 
holes is shown in Eig. Q for I up to 8 . Eor B = 0, the 
potential is Poschl-Teller, and becomes effectively wider 
and wider as B increases. 

The same step sizes are used as in section HI The 
results for greybody factors are shown for both small and 
near-extremal black holes in figures Q and ([^. 


V/k^ 



FIG. 5: Effective potential for kRh = O.OI and different 
values of I and coupling strength B. Two terms 
contributes, one due to the black hole geometry, i.e. Eq. 
© (left crests, in black), and the other due to coupling 
with background fields, i.e. Eq. (40) (right crests, in 
red). 


It turns out that effect of the coupling can be measured 

by a parameter B ^ | f (;£)', which 

is independent of the graviton mass. Note that a and kq 
are dimensionless parameters in massive gravity. In the 
small black hole case, 7 / is moved to higher frequencies as 
B gets larger due to higher potential barrier, similar to 
increasing 1. This argument also holds for near-extremal 
black holes. Besides, nonzero coupling will also make 7 ; 
show additional peaks and dips at some frequencies, for 
relatively large values of B and smaller values of I, for 
both small and near-extremal black holes. 

To explain the appearance of the visible peaks (and 
associated dips) in the near-extremal case, we write ex¬ 
plicitly the form of the effective potential appearing in 


^ A stepsize of 10 ^ is also sufficient for the numerical calculation 
for Z up to 8 and B up to 80, in the near-extremal case. 
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FIG. 6: Greybody factors for small black holes kRh = 0.01, with coupling between test scalars and background 
vector modes. The spikes in the greybody factors are due to resonances in transmission rate which in principle 
increase the transmission and allow more particles to penetrate the barrier. 


equation Q to the leading order in surface gravity. 


V{r) 


- 1(1 + 1 )- ^ - 

3 {kRh)'^ cosh^{kHr*) 

+0{kH^). 


k'jj{kHr*Y 

cosh^(fc//r*) 

(41) 


As B gets larger while keeping I fixed, the second term 
makes this potential barrier effectively wider, more sim¬ 
ilar in shape to the rectangular potential barrier or the 
double delta-function barrier(see Fig. Q), where the res¬ 
onances in transmission rate are well known [34] . A reso¬ 
nance in transmission rate increases the transmission, so 
more particles will penetrate the barrier, which in turn 
produces a spike in the greybody factor. This qualita¬ 
tively explains the resonance in the present scattering 
problem, as well as the tendency that the value of B 
needed for resonance grows with 1. In general, what kind 
of analytical potentials do or do not display the resonance 
in scattering is still an unsolved problem. 

Although the analytical form of the effective potential 
V{r) is not available for the case of small black holes. 


peaks are also expected when B is large enough, from 
the above discussion. This is verified by Fig. §. 


For better comparison with the blackbody radiation 
and the case without coupling, we plot the particle emis¬ 
sion rate {n(u})) and power spectrum dE/dtd\, as shown 
in figures ([To|)([l^. Although there are additional 

peaks in the greybody factors, which can raise the num¬ 
ber and power emission, these peaks moves to higher fre¬ 
quencies as the coupling becomes stronger and so become 
less important. For example, for small black holes with 
Th « 8k in Fig. ([^ for the B = 1 case, only the peak 
in 7o(a;) show up in figures Q(10l, but the peaks in 
7i(a;),72 (w), 73(0;) are all suppressed by the Boltzmann 
factor. 


Similar things happen in the near-extremal case, where 
Th = Spikes move to higher frequencies as B gets 

larger, reducing their effects on the emission rate and 
spectrum. However, as seen from Fig. there are 

also ’’incidental” spikes. For example, I from 0 to 8 all 
contribute spikes near oj/kn = 4. So, conservatively, we 
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FIG. 7: Effective potentials for near-extremal black holes, with various coupling strength between test scalars and 
the background vector fields, and I up to 8. Different couplings are differentiated by different colors. 


only calculate the emission rate and spectrum for B — 
0,1, 5,10, using the greybody factors from I up to 8. 

VI. CONCLUSIONS 

In this paper, we studied the greybody factor correc¬ 
tions to the Hawking radiation of black holes in massive 
gravity theory. We did the calculations numerically for 
small black holes, and analytically for the near-extremal 
ones (where the event and cosmological horizons are of 
the same order). Similarly, features of 'Yi{uj) were found, 
with very small values at low frequencies while approach¬ 
ing unity at higher frequencies, and the main contribu¬ 
tion to the radiated power comes from low-Z modes. 

To the best of our knowledge, the hrst results on grey¬ 
body factors of Schwarzschild-dS black holes (though in 
a different underlying context) were obtained by Kanti 
et al. [53]. Also, in recent papers ini isg, analyti¬ 
cal form of the greybody factors for Schwarzschild-dS 
black holes were obtained, consistent with our results, 
by matching the solution forms (usually hypergeometric 
functions) in different regions of r. However, these ap¬ 
proximations only hold either for small black holes, or 


for low energy modes. In contrast, we worked out the 
calculations for small black holes numerically, while for 
the near-extremal black holes analytically (exactly). This 
special case for the near-extremal black holes was found 
more than 30 years ago for normal mode calculations m, 
but remained unnoticed for computing the transmission 
coefficients. 


Once we couple the test field to the the background 
Stiickelberg fields, the situation becomes more compli¬ 
cated. The background fields can be decomposed into 
two modes, of which the vector mode has r-dependence 
only, while the scalar mode has t-dependence. We con¬ 
sistently couple the test field with the vector mode only, 
which in turn drives the greybody factors to higher fre¬ 
quencies and introduces some non-trivial features (reso¬ 
nances) as explained in Section [vj 

Recently, a solution for charged black holes in nonlin¬ 
ear massive gravity was found in |36j . It would be inter¬ 
esting to find the greybody factors for such black holes 
and compare them with those found here. 















7o 


7i 


72 




FIG. 8: Greybody factors for near-extremal black holes, with various coupling strength between test scalars and 
background vector modes. The black dashed lines are the analytical approximation of Eq. (291 for zero coupling, 
which can be seen to coincide with our numerical calculations. Different couplings are differentiated by different 

colors. 
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FIG. 9: The particle emission rate of small 
Schwarzschild-dS black holes {kRh = 0.01) with 
different coupling with the background vector fields. 

The amplitudes of the black hole emission have all been 
amplified by 10"^ times to be more visible. Different 
couplings are differentiated by different colors. Note 
that the greybody spectrum is always below the 
blackbody one before the artificial amplification. 

dE/dtdX 


300 



FIG. 10: The emission power spectrum of small 
Schwarzschild-dS black holes {kRh = 0.01) with 
different coupling with the background vector fields. 
The amplitudes of the black hole emission have all been 
amplified by 10 times to be more visible. Different 
couplings are differentiated by different colors. Note 
that the greybody spectrum is always below the 
blackbody one before the artificial amplification. 
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FIG. 11: The particle emission rate of near-extremal 
Schwarzschild-dS black holes with different coupling 
with the background vector fields. Different couplings 
are differentiated by different colors. The emission rates 
for 5 = 1,5,10 are amplihed by 500, 5000, 50000 times 
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